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Abstract 

A soft-wall warped extra dimension allows one to relax the tight constraints 
imposed by electroweak data in conventional Randall-Sundrum models. We 
investigate a setup, where the lepton flavour structure of the Standard Model 
is realised by split fermion locations. Bulk fermions with general locations are 
not analytically tractable in a soft-wall background, so we follow a numerical 
approach to perform the Kaluza-Klein reduction. Lepton flavour violation is 
induced by the exchange of Kaluza-Klein gauge bosons. We find that rates 
for processes such as muon-electron conversion are significantly reduced com- 
pared to hard-wall models, allowing for a Kaluza-Klein scale as low as 2 TeV. 
Accommodating small neutrino masses forces one to introduce a large hier- 
archy of scales into the model, making pressing the question of a suitable 
stabilisation mechanism. 
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1 Introduction 



Over the last ten years there has been a large increase in the study of extra di- 
mensional models following the realisation that they could help explain some of the 
unresolved problems in the Standard Model (SM). In 1999, Randall and Sundrum 
showed that a warped extra dimension could offer a geometric solution to the gauge 
hierarchy problem In the original Randall-Sundrum (RS) model, the fifth di- 
mension consists of a slice of AdS space bounded by ultraviolet (UV) and infrared 
(IR) branes. The warped space produces an exponential difference in energy scales 
between the two branes which solves the hierarchy problem. Matter fields were 
originally confined to the IR brane, however, it was soon realised that by allowing 
fermions to propagate in the extra dimension, the SM fermion mass hierarchy can 
be explained. By varying the location of the fermion wavefunctions in the fifth di- 
mension, the full scale of fermion masses from neutrinos to the top quark can be 
generated using only order unity parameters This setup also contains a 

built in mechanism suppressing unobserved flavour changing processes that result 
from couplings between SM fermions and excited gauge bosons which appear in the 
model 

Further interest in warped extra dimensions was generated by the AdS/CFT 
conjecture, when it was realised that the RS scenario is holographically dual to 
strongly coupled 4D field theories , Isf . It was in this context, studying AdS / QCD 
models, that the idea of a soft wall was first introduced \ld\. The soft wall is realised 
by removing the IR brane so the extra dimension extends to infinity, and replacing 
it with a smoothly varying spacetime cut off. The original AdS/QCD motivation for 
this was to more faithfully reproduce the linear Regge-like mass squared spectrum 
of excited mesons as opposed to the usual quadratic spectrum found in hard wall 
RS models. 

Inspired by the possibility of qualitatively different phenomenology, the soft wall 



scenario was subsequently applied to modelling electroweak physics ll2| . These 



models successfully showed that a soft wall extra dimension is generally less con- 
strained by electroweak precision observables than its hard wall counterpart, typi- 
cally allowing Kaluza Klein (KK) modes with masses of a few TeV. An important 
issue is related to the stability of a soft-wall setup, which is an open question in the 
models discussed in Refs. ll2|. Such a mechanism was suggested in Ref. 13 



promising the soft-wall extra dimension to equally well resolve the gauge hierarchy 
problem. 

With the removal of the hard-wall brane the Standard Model matter fields must 
necessarily propagate in the bulk. Graviton fluctuations and gauge fields were suc- 
cessfully analysed in this background but it was found that fermions presented par- 
ticular technical difficulties and only a simplified single generational model was de- 
veloped. Later studies of fermions in a soft-wall extra dimension have developed 



solutions to the fermion problem [IJ, |15|, |16| and have considered the experimental 



constraints imposed by the electroweak observables. However, the fermion flavour 



1 



pattern of the SM has not been considered in much detail, in particular with respect 
to the generation of neutrino masses and the experimental bounds on lepton flavour 
violation. 

In this paper we present a numerical solution to analyse a single generation of 
fermions in the soft-wall extra dimension. We extend this solution to three genera- 
tions by treating flavour mixings as perturbations to the original solutions, and apply 
it to the lepton sector of the SM. We construct a setup, where the lepton flavour 
pattern is accommodated by flavour dependent localisations. It is shown that in 
order generate small Dirac neutrino masses by this mechanism we need to introduce 
a hierarchy of scales of order 10^^ into the model, making crucial the issue of a 
suitable stabilisation mechanism. We finally carry out an analysis of the constraints 
coming from various lepton flavour violating processes, averaging over random order 
unity Yukawa couplings, and find that models with only a modest hierarchy of scales 
are relatively mildly constrained, whereas the model with a large hierarchy allowing 
sub-eV neutrino masses lies well within current experimental constraints, even for 
a KK sea 1^ of 2 TeV. In the latter, flavour violation is considerably suppressed 
relative to its hard wall counterparts, such as the ones analysed in js], [l^ , and the 
range of masses lies in the reach of the LHC experiment. 

At this stage we do not try to accommodate the flavour structure of the quark 
sector, which should be possible in a similar way. Also we reproduce the neutrino 
masses and mixings only at the qualitative level, which is sufficient to estimate the 
rates of lepton flavour violation. 



2 Bulk Fields in a Soft- Wall Extra Dimension 



Our conventions follow most closely those laid down in Refs. [12|, |16|. The 5D 
spacetime has metric 

els'" = e-^^^y^MNdx'^dx'', (1) 

where y represents the extra spatial dimension and timn = diag(+, — , — , — , — ). We 
take a pure AdS metric, A{y) = log ky, where k is the AdS curvature scale. There is 
no IR brane, the extra dimension extends to infinity and the soft wall is introduced 
via a dilaton field $ with the action describing gauge and matter fields given by 

//•oo 
d'^x / dy^e-^C. (2) 
Jyo 

Here y^ = 1/k is the location of the UV brane. The dilaton field is taken to have 
the following power law behaviour 

Hy) = M'- (3) 



^When referring to the soft wall model we take the KK scale to be the mass of the first KK 
mode of the Z boson. 



2 



The dimensionful parameter /i will set the mass scale of the lightest KK excitations. 



The behaviour of other powers has been discussed in detail in Ref . . It is shown 



in Ref. 12[ that an appropriate form for the Higgs VEV in such a background is 
given by 

hiy) = r^e/'^iV (4) 
where 77 is a dimensionless 0{1) coefficient. 

2.1 Massive Gauge Fields 



In Refs. [12|, |16[ only massless gauge fields are considered as the gauge couplings 
being considered are assumed to be between quarks and massless gluons. The fiavour 
changing neutral currents (FCNC) we will be considering here are mediated by the 
massive Z boson and so we first develop the solutions for such a gauge field. A 
massive gauge field propagating in the bulk has the Lagrangian 

C = -\fmnF'''' + ^MIAmA''. 

The mass term for the weak gauge bosons arises from spontaneous symmetry 
breaking and with the Higgs VEV as given in (jll) we have = \g1h{iiY. Varying 
the action we obtain the equation of motion 

v g 

Imposing the gauge Ay = 0, inserting the Kaluza Klein (KK) reduction, 

00 

/l,(x,y) = 5^4")(x)/t)(y), (5) 

n=0 

and requiring the A^^\x) to be mass eigenstates, we find the have to satisfy 

and are canonically normalised by 

yo 

A complicated analytic solution of the above equation of motion was developed in 



14| but for our purposes we find it more convenient to solve the equation of motion 
numerically. We apply Neumann boundary conditions to the wavefunctions at yo 
and vary rjg^ in order to find a normalisable solution with the appropriate 4D zero 
mode mass of mo — 91 GeV for the Z boson. 
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It is interesting to note that unlike in hard-wall models, the profile of a massive 
gauge boson is independent of the curvature scale k. This can easily be seen by 
looking at the equation of motion ([6]) where k only appears in the term involving 
the 5D mass Ma which is the term that must be varied in order to generate the 
correct zero mode mass. The profile of the first few KK modes of the Z boson 
are plotted in Fig. 1. In Fig. 1(a) the zero mode and the first two KK modes are 



plotted with respect to a flat metric. Fig. 1(b) shows the proflle of the zero mode in 



the form that it couples to the fermions (see later). The UV behaviour of the zero 
mode is less flat than in hard wall models and could possibly lead to large violations 
of universality of gauge couplings once fermions reside at different locations in the 
extra dimension. Note also that as is the case for massless gauge bosons 12|, |l6| the 
higher KK modes become more and more IR localised, a fact that will be important 
when considering couplings between higher gauge modes and fermions. 

The KK spectrum for the Z boson for different values of /i with k = 10'' TeV is 
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Note that the mass of the flrst KK mode mi ~ 2/x, and hence the KK scale, 
Mkk, scales with /i. Higher modes follow a Regge hke spectrum ~ n. 
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Figure 1: Z boson proflles for k = 10^ TeV and = 1 TeV. (a) The zero mode 
(solid), flrst (dashed) and second (dotted) KK modes plotted with respect to a flat 
metric, (b) The zero mode proflle as it couples to fermions. 



2.2 Fermions 

We consider 5D Dirac spinors and which are components of doublets and 
singlets under SU{2)i respectively. Note that L and R do not denote chiralities. 
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but are related to the charges under SU{2)l. The chiral projections of these spinors 
are l± = |(1 =f 7^)\1//^, same for The action for two free fermions in the bulk 
is 



Jyo 



The fiinfbein and spin connection for the metric ([T]) are = €^^"^^5^ and ojm = 
(—^7^7^, O) and the covariant derivative is then Dm = Qm + ^m- M^ r are the 5D 
Dirac masses related to 

The difficulty of placing uncoupled fermions in the soft wall background is well 
documented in Ref. [l2|, all solutions suffer from divergent gauge couplings for high 
enough KK modes. The underlying problems stem from the non compact nature 
of the extra dimension and it is shown that in order to find workable normalisable 
solutions, the Yukawa couplings between fermions and the Higgs must be taken into 
account. An alternative approach to introducing Yukawa couplings was presented 
in Ref. 15|] where a y dependent Dirac mass term is introduced, somewhat like the 
y dependent bulk mass arising from the Higgs VEV in the case of the massive gauge 
boson above. Here, however we will stick to constant Dirac mass terms Ml and Mji 
and introduce Yukawa couplings into the action: 

//•oo ^ 
d^x / dyy/ge'"^^ [ ^^(x, y)h{y)-<^ r{x , y) + '^r{x, y)h{y)'<^ l{x, y 

Defining ^l,r = ^^'^^^^'^'4^l,r and m(y) = ^h{y), the equations of motion are 

il^d^il)L,R± ± dyijL,RT - e'^M^Rip^RT - e''^m{y)'ipR^LT = 0- 
Using the KK reduction 

oo 
n=0 

and requiring the il)^jf^{x) to be mass eigenstates, the /"^"^s will be given by 




Ml m{y) 
m{y) Mr 





(7) 



The V'j!^^ will be canonically normalised by 

/•oo 

/ dy (ft' m + ft2 ftl) = Smn. 
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We impose Dirichlet boundary conditions at uq for /|"^ and /jj*^ in order to obtain 
a chiral 4D theory. 

Analytic solutions of ([7]) are only possible for a small set of Dirac mass terms, 



namely Ml = Mr and Mi+Mjiiik = [16|. As in hard- wall models, the Dirac mass 



parameters Ml^^ dictate how the fermion is localised in the extra dimension and 
it is convenient to parametrise these in terms of the AdS curvature, M^^ji = CL^^k. 
Unfortunately the sets of parameters for which analytic solutions are available do not 
explore the full geography of possible mass parameters and as we shall see may lead 
to situations with unacceptably large rates of flavour violation. Ideally we would 
like to solve ([7]) for any set of Dirac masses and this requires a numerical approach. 

The numerical solution we have developed involves a shooting type method. In 
order for the solutions to be normalisable they must not diverge in the IR and this 
only occurs for the correct choice of m„, thus generating the KK spectrum. We 
choose a suitably large distance L into the IR and solve the equations of motion 
subject to the UV boundary conditions and a starting choice for m„. We then 
iterate the solution using Newton's method in order to find a value for m„ such that 
two of the solutions e.g. fl^_^ and converge to zero at y = L. The equations of 
motion then automatically ensure that the other two solutions will also converge to 
zero for large y. Our solution has the advantage that it seems to be quite capable 
of finding solutions even for large values of the AdS curvature scale k, however it is 
not so suited to finding solutions for multiple generations of fermions as is done in 



Ref 16 



3 Leptons 

3.1 General Considerations 

Due to the presence of the extra KK states in extra dimensional models, couplings 
between SM particles and their KK excitations can potentially lead to conflict with 
experimental observations. In the SM, a tight set of constraints comes from the 
experimental bounds on flavour changing neutral currents (FCNC). In the hard- 
wall Randall Sundrum model these processes have been investigated and are shown 
to occur at rates that are dependent on the fermion locations. However, there are 
certain choices of fermion locations which provide almost universal gauge couplings 
and these almost universal gauge couplings are the source of the so called RS GIM 



mechanism which suppresses FCNC [3| 



It was found in Ref. [12] that for fermions in the soft-wall background, the 
analytic solution with Ml = Mr can produce a large hierarchy of masses but only 
one of the fermion pair {ip^ , ) could reside in an area of universal gauge couplings 
and it was thus assumed that dangerous rates of FCNC would be generated in such 
a situation. In fact this is one of the main motivations for finding a numerical 
solution to the fermion equations of motion, in the hope that one would be able to 
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find fermion locations which can give a large mass hierarchy and yet simultaneously 
reside in an area of universal gauge couplings. 

The gauge interaction between bulk gauge bosons and fermions is given by 

//•oo 

The couplings of ip^^ to different KK gauge modes is then 



9l = 95 I dy 

yo 



2 ' (0) 



(8) 



The dependence of the gauge couplings on the fermion locations cl^r is shown 
in Fig. 2. It can be seen that for cl > 1/2 the couplings become universal. In the 
case where cl = cr, the couplings of one of the fermions would lie in the universal 
region the other would lie in the opposite part of the plot. However, with opposite 
Dirac masses, cl = —cr, we are able to place both fermions in a region of universal 
coupling at the same time and we would thus hope to suppress FCNC. 

In hard wall models the origin of the regions of universal couplings is quite clear 
and derives from the profile of the gauge field wavefunctions which are fiat in the 
UV. Hence if fermion profiles are relatively UV localised the gauge couplings will 
be universal. However, in the soft wall model, looking at the profile of the zero 



mode of the massive gauge boson in Fig. 1(b) it is certainly not flat and one may 



wonder why we still find regions of universal gauge couplings. The explanation can 
be seen by considering the fermion profiles. Fig. [3] shows the fermion wavefunctions 
contributing to the gauge coupling of ip^ for cl = —cr = 0.7 which are locations 
that live in an area of universal gauge couplings. Whilst f^^ is heavily UV localised 

we see that f^l is actually peaked into the IR which we would expect to contribute 
to non-universal couplings. However, when one considers the relative size of the 
contributions of each of these wavefunctions to the gauge coupling as given by Eq. 
[8] we find ~ 3 x 10~^ i.e. almost the entire contribution to the 

gauge coupling comes from f^^ which is heavily UV localised. Although the relative 

dominance of f^^ is not that clear to see from Fig. [3], it becomes obvious when one 
realises that it has a value at fiy^ of about 2000. Also, its extreme UV localisation 
can be seen by the fact that 99% of the area of (f^l)'^ lies in the region /xy < 0.01. 
Hence the dominant contribution to the gauge coupling comes from a region in the 
extreme UV where the gauge profile is effectively given by its UV boundary value, 
thus producing universal couplings for fermions. 

In order to generate the fermion mass hierarchy seen in the SM we have to 
carefully choose the c parameters. The zero mode masses for different c parameters 
can be seen in Fig. |H Unfortunately, the shape of the plot presents a problem for 
simultaneously generating a large hierarchy of masses and universal gauge couplings. 
Whilst it is easy to generate a large hierarchy of masses for the choice of parameters 
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(a) (b) 

Figure 2: Gauge couplings of the Z boson (a) and its first two KK states (b) as a 
function of the fermion location, normalised so that the coupling to the zero mode 
is unity for large cl- k = 10^ TeV, /i = 1 TeV and cl = —cr. 

cl = cr, as has been stated above, this is likely to lead to high rates of FCNC. In 
order to avoid these unacceptable rates we would like both the fermions to reside 
in an area of universal gauge couplings, this corresponds to the top left corner of 
the contour plot where we have cl > 1/2 and c/j < — 1/2. In this area the zero 
mode mass bottoms out at around fi'^/k and it is not possible to create a large 
mass hierarchy. The solution to this is to simply increase the hierarchy of scales in 
the model. Keeping /i = 1 TeV we can see from Fig. [5] that with k/n = 10'' and 
Cl = — c/j the zero mode masses could cover the full range of charged lepton masses 
whilst remaining in an area of universal gauge couplings. 

The reason for the zero mode mass having a minimum value for cl = —cr can be 
seen by considering the two different ways small masses are generated in such models. 
The zero mode masses are generated via Yukawa couplings which involve the overlap 
between ^ r and the Higgs VEV. In the case where = cr, the wavefunctions 
of the zero mode fermions become oppositely localised and can be arranged to have 
an arbitrarily small overlap with each other thus creating arbitrarily small masses. 



this is the mechanism used in the "split fermion" model [18| . However when we take 
Cl = —cr the fermions completely overlap each other and the zero mode mass is 
then entirely determined by their overlap with the Higgs. In RS models where the 
Higgs resides on the IR brane, this overlap can be made arbitrarily small by heavily 
localising both the fermions at the UV brane. However in the soft-wall model where 
the Higgs must necessarily propagate in the bulk and has a non zero value at the 
UV brane, the fermion wavefunctions will always have a minimum overlap with the 
Higgs dX yQ = l/k however much they are UV localised. 
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Figure 3: Fermion zero modes with cl = —cr = 0.7 for k = 10^ TeV and /i = 1 



TeV. solid, f^^ dashed. Note that f^_l takes a value of around 2000 at yo 
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Figure 4: Contour plot of log^g (^o//^) for the zero mode masses of fermions with 
/X = 1 TeV and k = 10^ TeV. 
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Figure 5: Fermion zero mode masses with cl = —Cr, /i = 1 TeV and k = 10^ TeV. 



3.2 Neutrino Masses 

We can generate small Dirac neutrino masses by introducing right-handed neutrinos 
and allowing Yukawa couplings between them, the Higgs and the left-handed neu- 
trinos. A similar approach has been taken in the hard-wall case in Refs. gEi- The 
left-handed neutrinos will share the same cl parameters as the corresponding left- 
handed charged leptons since they are part of the same doublet under SU{2)l. We 
are then free to place the right-handed neutrinos in a suitable location in order to 
generate sub-eV masses. However, since we still require cl > 1/2 we again find that 
we are unable to generate such small masses without vastly increasing the overall 
hierarchy in the model. It seems necessary therefore to work with a hierarchy similar 
to that proposed in the original Randall Sundrum model. The issue of stabilising 
such a large hierarchy is an important question, and it would be very interesting to 
redo our analysis in context of the stabilised model proposed in Ref. ISj. With this 
in mind we choose k/jj, = 10^^ and are able to produce neutrino masses of order 0.1 
eV by choosing cl = 0.6 and cr = —1.3. 



3.3 Three Generations 

When incorporating all three generations of leptons into our model the Dirac mass 
terms Ml and Mr and the Yukawa coupling constants are promoted to 3 x 3 matrices 
mixing the different generations. We assume that the basis of states in which the 
Ml and Mr are diagonal does not correspond to one in which the Yukawa couplings 
are diagonal. Rather than finding exact solutions for all three generations in such 
a scenario, our approach to this problem follows closely the method used in Ref. jsj 
for the Randall Sundrum model. We solve the equations of motion individually for 
each generation with a Yukawa coupling A5 = 1, excluding fermion mixing, and use 
these solutions as basis from which we treat the full matrix of Yukawa couplings 
including mixings between the generations as perturbations. We specifically choose 
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a large number of random Yukawa couplings, taking ^ < \X5ij\ < 2 with random 
sigiicl, and require that the average zero mode masses reproduce the observed lepton 
masses. We also choose to locate the left-handed fields of each generation close to 



each other in order to generate large neutrino mixings, in the spirit of Ref. |19 
However, we do not aim at reproducing the neutrino masses and mixings precisely. 
All we arrange for is an overall neutrino mass scale of order 0.1 eV. Using our tools, 
a full model of neutrino masses could be constructed. However, for the following 
estimate of lepton flavour violation these details are not needed. Also we use the 
fact that neutrino mixings are order unity. 

In the case where we are not interested in generating neutrino masses via loca- 
tions, we take only a moderate hierarchy of scales, k/fi = 10"^. In this regime we 
choose the following three scenarios: 

(A): CLi 



(B) : CLi 

Cm 

(C) : CLi 



In the regime where we can also generate neutrino masses, k/ jj, = 10 we choose: 

(D): CLI = 0.60, CL2 = 0.60, cl3 = 0.60, 
cm = -0.82, c^j,2 = -0.64, = -0.55. 

Our choices for the different scenarios (A), (B) and (C) are to demonstrate the 
effects of degenerate c^ localisation (A), small separation in the c^ to introduce some 
non-universality in the left handed sector (B), and placing the left handed fermions 
closer to the IR brane (C). We expect the behaviour to be quite general and thus 
only choose one scenario with a larger hierarchy. The mass of the first fermion KK 
states is about 1.5 TeV. 
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= -0.790. 



4 Flavour Violation 

With a full three generations of leptons implemented as above, the transformation 
to fermion mass eigenstates will induce flavour violating couplings to gauge fields, in 
particular the Z boson and its KK excitation^ We define the neutral current gauge 

^We do not consider CP violation, i.e. we take As^ to be real. 

^Note that we work in a gauge field basis, where the Z boson zero mode is massive. Its properly 
weighted wave function is not flat, resulting in non-universal couplings to fermions. Alternatively, 



11 



couplings in the basis of mass eigenstates as 

where the unitary matrices IA± diagonahse the full fermion mass matrices and Q^'^ 
are diagonal matrices that contain the couplings of the nth KK state of the Z boson 
to each fermion state as derived from Eq. [8] and normalised to the coupling of the 
muon (see also Ref. jl]). Flavour violation induced by these couplings is dependent 
on the non-universality in the couplings of different flavour states and the mixing 
between the states. Different fermion locations increase the non-universality but at 
the same time lead to small mixing angles. Conversely similar fermion locations 
produce large mixing but this is compensated by universal couplings. 

As was done in Ref. |5|] we calculate the rates of the various flavour violating 



processes using the techniques developed for family non-universal Z' bosons [20 
The main difference being that there is no mixing between the different KK states 
of the Z boson, while the zero mode also has flavour violating couplings. 

The first process we consider is the tree level exchange of a Z boson and its KK 



states mediating the process Ij — )■ liUli. The rate for this process is given by [20 

dm' 

t]\ 'I tj\ 'I ij 



where 



^5 = E^(^i"^)^.-(^?^) 



M2 

n " 

where we take the sum over the Z boson zero mode and the first two KK modes. M„ 
is the mass of the nth KK mode of the Z boson. Due to the Regge type behaviour 
of the KK spectrum it is not clear that the above series should converge. However 



as noted in Ref. [16|, due to the increasing IR localisation of higher gauge boson 
KK modes, the couplings rapidly decrease and the series converges after only a few 
terms. In the fermion sector we take into account only the zero modes. Mixing with 
KK fermions is small, leading to negligible effects at the current precision. 

The branching ratios for the above processes in the different scenarios we consider 
are then found to be 



Br(yU — )■ eee) : 
Br(r — )■ yUyU/i) : 
Br(r — )■ eee) : 



(A) 




(B) 




(C) 




(D) 




2.7 X 10" 


-14 


5.1 X 10- 


-12 


4.6 X 10- 
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2.4 X 10" 
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7.0 X 10" 
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2.6 X 10" 
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-13 


2.8 X 10- 


-16 



one could work in a basis, where the zero mode is massless. Then only the KK states would couple 
non- universally. 
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These numbers are obtained for a KK scale of 2 TeV. They are the result of averaging 
over random Yukawa couplings in the range stated above. The experimental bound 



Br(/i — )■ eee) < 1.0 x 10~^^ 2l|] is satisfied in the cases (A) and (D). However, it 
appears that the couplings are not universal enough to allow for much separation 
between the left handed states (B), and placing the fermions too close to the IR 
brane (C) also exceeds the experimental bound. However, like in the hard-wall case, 
the rate for this process depends on the KK scale as 1/M|^^. Thus with a KK scale 
of twice as big (i.e. take fi = 2 TeV while keeping k = 10^ TeV) scenarios (B) and 
(C) would also acquire an acceptable rate. The experimental bounds for the other 
two processes Br(r — )■ /x/i/2) < 2.1 x 10~^ and Br(r — )■ eee) < 2.7 x 10~^ Q are 
well satisfied in all the scenarios. Note that in Ref. [sl in a case similar to (A), (C) 
and (D) a branching ratio Br(yU — t- eee) = 5 x 10~^^ has been found for a KK scale 
of 10 TeV, translating into Br(yU — )■ eee) = 3 x 10~^^ for a KK scale of 2 TeV. This 
demonstrates that lepton flavour violation is suppressed by up to four additional 
orders of magnitude in the soft- wall case. 

We are also able to calculate the expected rate of /i — )■ e conversion in a muonic 
atom. The most stringent bound comes from the Sindrum-H Collaboration [23] in 
22Ti where Br(/i~iV — )■ e~N) < 6.1 x 10~^^. We can calculate the branching ratio 



for this process by 2ol 



^CAPT 



Br(^-A^ ^ e'N) = ' ^ ^ -f\Fp\' ( 1^-^ + 



where 



n " 

and we take 

Kln=9^L.., =0.19^7^,,, and B^^l^ = O.U g^^^^. 

Here, g^jj^^ ^ are the usual Standard Model quark couplings and we have taken the 
approximate values of the quark couplings to the higher KK gauge modes from the 



values derived in Ref. [16[. We also take Zeg = 17.6, Fp = 0.54 and Fcapt = 
2.59 X 10^ s"^ |2J], where Zeff is an effective atomic charge, Fp is a nuclear matrix 
element and Fcapt is the muon capture rate. 

We find branching ratios for the different scenarios of 



(A) (B) (C) (D) 

BT{fi-N -> e-N) : 1.6 x 10^^^ 3 3 ^ iQ-n 2.7 x lO^^^ 1.4 x 10"^^ 

Again we find scenarios (A) and (D) lie within the experimental bounds but 
separating the states (B) or placing them too close to the IR brane (C) produces 
an unacceptable rate. Again we find flavour violation suppressed with respect to 
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the hard- wall case [5|, making a KK scale of 2 TeV consistent with observations. 
However, next-generation experiments, such as PRISM at JPARC with a reach of 
Bt{ij-N e'N) ~ 10"^^ - 10"^^ could probe a KK scale of 6 - 20 TeV, i.e. the 
interesting parameter range of the present model. 

A third set of processes considered in Ref. 2^ are one-loop radiative lepton 
decays. Here the decay width is 

aG%mf /, ..,2 , ..,2\ 

ni.-^ia) = ^^{\eif + \e\'), 

where the dipole moment couplings of an on-shell photon to the chiral lepton currents 
are given by 



«2 = E]ilK"'"'^ei"' 



5 



where mi is the charged lepton mass matrix. Using this we obtain the following 

rates for radiative decays 

(A) (B) (C) (D) 

Br(^ 67) : 2.0 x lO^^^ 3.7 x lO^^^ 3.4 x lO^^^ 1.2 x 10"^^ 
Br(r -> /i7) : 5.2 x lO^^^ 1.9 x 10-^^ 9.4 x lO^^^ 4.1 x lO'^^ 
Br(r ^ 67) : 3.7 x 10"^^ 1.5 x 10"^^ 9.2 x lO^^^ 2.6 x 10"^^ 

All of these branching ratios lie well within the experimental bounds Br(/x 
67) = 1.2 X 10-^1 |25j, Br(r ^ /i7) = 4.4 x 10"^ and Br(r ^ 67) = 3.3 x 10"^ 



Again these rates are suppressed relative to their hard-wall counterparts. 

In scenario (D), processes such as — )■ 67 can also be mediated by the KK states 
of the sterile neutrinos. This process was investigated for the Randall-Sundrum 
model in Ref. ji^]. We use the formalism developed there and we find that in our 
model the branching ratio for this process is given by the relative coupling strength 
of the muon, the W and the KK muon neutrino to the zero mode muon neutrino 
times a loop factor. Because of the large mass differences of the sterile KK neutrinos, 
the GIM mechanism breaks down. We assume that the neutrino mixing angles are 
large, not leading to any suppression of the rate. Given a relative coupling of 0.0057 
we find Br(/i — )■ e7) = 1.5 x 10~^^ which again lies within the experimental bounds 
for a KK scale of 2 TeV. 

There are also contributions to /i — )■ 67 related to the exchange of KK fermions 



28L |29| . which were neither included in our estimate above nor in Ref. [27j. These 
contributions dominate the rate of radiative lepton decays in hard-wall models. A 
similar behaviour is likely in the soft-wall model. However, given the suppressed 
rates for flavour violation in the latter, we expect that even including these extra 
contributions, the rate for /z — )• 67 and a KK scale of 2 TeV should not exceed the 
experimental bound. 

We have shown that lepton flavour violation can be suppressed in the soft wall 
model and the amount of suppression depends on the fermion locations and the 
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hierarchy of scales, k/ ^. The reason for this suppression comes from the abihty to 
place the fermions in regions of universal gauge couplings i.e. heavily UV localised 
as explained in Section 13.11 When the hierarchy of scales in the model is increased 
by fixing and increasing k the gauge boson profiles remain unchanged whilst the 
fermion profiles become even more UV localised thus providing an even greater 
suppression of flavour violation. 

5 Conclusions 

In this paper we have studied the lepton sector of the SM in a soft-wall extra di- 
mension, applying flavour dependent fermion locations to accommodate the observed 
lepton flavour structure. The Higgs is a bulk field, with a VEV that increases near 
the soft wall. We have in particular considered the inclusion of small Dirac neutrino 
masses and investigated the constraints on the model from lepton flavour violation 
mediated by the Z boson and its KK states. In order to do so we first developed so- 
lutions for a massive gauge boson in the soft-wall background and found the profile 
is independent of the AdS curvature scale. In order to generate the masses of the 
charged leptons whilst keeping the fermions located in an area of almost universal 
gauge couplings we find that we need to increase the hierarchy of scales in the model 
to around k/fi = lO''. When incorporating sub-eV neutrino masses we need a much 
larger hierarchy, and we choose k/fi = 10^^, similar to the hierarchy between the 
Planck and the electroweak scales. 

To incorporate three generations of leptons into our model we solve the fermion 
equations of motion numerically, including an order one flavour diagonal Yukawa 
coupling and use these solutions as a basis of states from which we treat off-diagonal 
Yukawa couplings, connecting different generations, as perturbations. The mass 
term related to the diagonal Yukawa coupling is necessary to generate a normalisable 
wavefunction and cannot be treated as a perturbation. We can construct the full 
lepton mass matrices, including KK states and diagonalise them to find the fermion 
masses and mixings. However, to our level of precision we can neglect the fermionic 
KK states. The locations of the left-handed fermions are dictated by the fact that 
we require large mixings in the neutrino sector. We take a large number of random 
Yukawa couplings and choose the locations of the right-handed fermions so that the 
averaged zero mode masses reproduce the SM charged lepton masses. 

With the inclusion of off-diagonal Yukawa couplings, the transformation to mass 
eigenstates produces flavour violating couplings. We calculated the expected rates 
for various flavour changing processes for a number of different scenarios. We found 
that the soft-wall model is in fact mildly constrained when we consider a scenario 
with a low hierarchy of scales such as k/fi = 10^. The most stringent constraint 
comes from — )■ e conversion in a muonic atom where we find that only the scenario 
where all the left-handed leptons have degenerate locations well toward the UV brane 
would occur at acceptable rates with a KK scale of 2 TeV. This is a considerable 
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suppression of lepton flavour violation compared to hard-wall models, such as the 
one studied in Ref. jsf. Including a larger hierarchy of scales {k/ jj, = 10^^), it is also 
possible to generate sub-eV Dirac neutrino masses. In this case the model is even 
less constrained and most of the FCNC processes would occur at rates well below 
the experimental bounds. The most stringent bounds are coming from radiative 
decays, such as — )■ 67. Again a KK scale of 2 TeV seems sufficient to keep the 
rate below the experimental bound. Our estimate for this rate does not include 
contributions from KK gauge bosons, and it would be interesting to include these in 
a more detailed analysis. Another obvious direction of research would be to extend 
the present setup to the quark sector, similar to an a naly sis that was performed 



recently in much detail for the hard-wall model in Ref. [30 . 

The soft-wall extra dimension continues to offer a valid model for electroweak 
physics, with constraints from precision data relaxed compared to the hard-wall 
model. Having said this, we have found that with a (gauge boson) KK scale of 
2 TeV the complete lepton flavour structure can be accommodated while keeping 
rare processes below experimental bounds. In our setup the KK states of fermions 
have masses around 1.5 TeV, within reach of the LHC experiment. Thus the soft- 
wall framework seems to offer an alternative when it comes to suppressing flavour 



violation to models relying on flavour symmetries [29|, |3ll, |32| , a bulk Higgs |33i] or to 
utilising non-minimal representations under the SU{2)r bulk gauge symmetry [s^]. 

The parameter range with a large hierarchy k/^ = 10^^ is both attractive to 
further suppress flavour violation and necessary to accommodate neutrino masses. 
This rises the important question whether such a hierarchy can be stabilised, like in 



the way proposed in Ref. [13||. It would be very interesting to extend our analysis 
to such a framework. 
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